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Abstract — This paper present mathematical modelof local and global integral for elastic plastic damage 
mechanism of displacement discontinuity. The discrete nature of stress intensity factor approach where there is 
large scale yielding is inadmissible; hence, continuous energy path contour analysis becomes crucial. The 
simple energetic path can easily be evaluated computationally and the most prevalent method. However, fails to 
account for dynamic effect, kinetic energy and thermal effect, body force and effect of surface traction. The 
research is focused on deformation contour analysis of elastic-plastic void growth. In this paper, local 
andglobal representation of deformation field near a crack tip due to loading was examined using Divergence 
theorem and energy balanced method. The local integral is equivalent to energy release rate for nonlinear 
elastic materialunder quasi-static condition. The global J- Integral captured thermal effect, body force non zero 
surface traction and effect of plasticity. The integral permits more accurate analysis possible for all 
deformation fields including very close proximity of crack tip and all categories of material. TheFatigue crack 
growth analysis was developed to support economic fail safe and damage tolerance philosophies, and the 
execution, effectuation of damage tolerance control procedures for the workings of structural systems. 
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I. INTRODUCTION 

Plane strain deformation field surrounding the tip of a crack 
is of immense significance in damage mechanism of solid 
component in structural design and analysis. It deals the 
studies of the behaviour of stress and strain field in the 
neighbourhood of fissure tip. The product of stress and strain 
exhibits a singularity which is inversely proportional to the 
magnitude of crack from the tip and near the tip where there 
is numerous invisible deformation fields or contours. This 
represents energy line integral with the characteristics of path 
that does not depend on all the contours encircling fissure tip. 
In Rice (1968) this line integral has same value for all paths 
around the tip of a crack in two dimensional deformation 
field of linear and non-linear elastic material. The correct 
choice of this contour is crucial; the result gives approximate 
estimate of the degree of strain concentration near the tip of a 
crack in a notch crack problem. Irwin (2004) posited that the 
principle of continuous deformation path is related to path 
independence if a path cn was obtained from c 1 by 


continuously moving c 1 with ends fixed until it coincide 
with c 2 • This revealed thatthere are two of the indefinite 
numerous intermediate paths for which the integral always 
retain its value. Besides, there is a continuous deformation 
path of an integral keeping the ends fixed since deformation 
path always consist of only points at which the function is 
analytic, hence the integral retain same value. 

Contour integral has great contribution in fracture 
mechanics; it characterizes the parameter used for nonlinear 
materials by idealizing the elastic plastic deformation as 
nonlinear elastic.According to Boulenouar et al (2013) the 
use of crack propagation laws based on stress intensity 
factors is quite good in scientific and engineering application 
of mechanics of fracture. The stress intensity factor 
sufficiently defines the stress field close to the crack tip and 
provides fundamental information on how the crack is going 
to advance. Furthermore, Rice (1967) applied deformation 
plasticity to analyze crack in a nonlinear material; therefore, 
in nonlinear energy released rate, J integral could be used 
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and represented as a path independent line integral.Riceand 
Rosengren (1968) also showed that J integral uniquely 
characterize crack tip stress and strain in nonlinear materials. 
So.the J integral could be viewed as both an energy 
parameter and a stress intensity parameter. It is again seen as 
path independence contour integral for the analysis of crack 
problems. The contour integral could be evaluated along a 
contour surrounding the crack tip where numerical accuracy 
is better. In path independence, it is possible to compute J at 
remote contour provided an appropriate correction term that 
is an area integral is applied. Furthermore, for J integral to 

1 

remain path independent, the stress strain must vary as — 

r 

near crack tip. This energy release rate measuresthe energy 
available for an increment of crack extension also known as 
crack driving force or crack extension force. It has the 
potential ability to evaluate the integrity of structural 
component: it has become common to evaluate the integrity 
of structures using path independent J- integral interpreted as 
the intensity of the elastic - plastic deformation and the 
stress field surrounding a crack tip. However, the integral 
loses the property of path independence and can no longer be 
interpreted as the energy released rate under some 
conditions. It is not evident that if the elastic plastic stress or 
strain fields near a crack tip is characterizable in terms of J - 
integral. In addition, Oyesanya (2007) found that the integral 
cannot characterize the crack tip stress strain analysis of 
structures under thermal loading outside HRR dominance. In 
small size dependent fracture analysis characterized by strain 
gradient elasticity the convenient strain gradient J- integral is 
incapable of accounting for the deformation energy 
distribution from the strain gradient. It underestimates the 
critical load and over estimates the critical deflection then 
higher order J- integral was necessary for analysis of fracture 
behaviour of structures with none negotiable strain gradient 
Lam e’tal (2004). 

Due to complexities in investigating crack tip stress and 
strain of displacement amplitude loading condition, the 
discrete nature of stress intensity factor approach becomes 
inadmissible where there is large scale yielding.The need for 
continuous energy path contour becomes crucial. However, 
the simple energetic path fails to account for dynamic effect, 
kinetic energy, thermal effect and body force. It loses its 
physical significance as a crack driving force when 
deformation is not reversible. Therefore, in this paper, a 
global representation of deformation field near a crack tip 


due to loading was examined using globalized energy 
balance approach. It is a general representation of J integral 
which includes the existence of a fracture process region and 
the effect of plastic deformation, body force, thermal strain 
and inertial effect of material.Crack tip singularities are 
investigated using energy line integral that exhibits path 
independence for all contours near the tip of the crack in a 
two dimensional deformation field of a linear elastic 
material. Besides, the global J- Integral is used to extract the 
magnitude of crack tip stress intensity factor in fatigue crack 
problem and the characteristics of the J-integral under mode- 
1 loading. The integral permits accurate analysis possible for 
all deformation fields including very close proximity of 
crack tip. This arbitrary continuous and differentiable 
contour represent the magnitude of energy available for 
crevice propagation, that is, the value of this integral equates 
the energy released rate of linear and nonlinear elastic plastic 
body which contains void. 

II. LOCAL ENERGETIC PATH CONTOUR 
INTEGRAL 

Plastic zone size plays crucial role in actual crack 
propagation, hence elastic fracture dynamics that can be used 
notably among is J-integral. It is a parameter that deals with 
energy while stress intensity factor is a local parameter that 
deals with displacement and stress field in the vicinity of a 
crack. The approaches may be different but the goal is the 
same to characterize a crack. The magnitude of plastic zone 
affects that actual fatigue crack progression. Therefore, 
elastic plastic fracture dynamic model for crack growth 
becomes crucial. This gives rise to energy flux integral that 
denotes path independent energy line integral that quantifies 
the strength of the singular stress and strain field in the 
vicinity of a crack tip (a measure of energy available for an 
increment of crack extension, the energy released rate). It is 
also known as crack extension force or crack driving force. 
The magnitude of this J-contour integral can be calculated 
from a path independent contour integral in the 
neighbourhood of the flaw tip. 

However, taking into view an arbitrary contour cthat 
encloses a crack tip whereby the path is in anticlockwise 
direction to evaluating J along the contour, the value of this 
integral called contour integral. This is equal to the energy 
released rate in a nonlinear elastic material that contains a 
crack.In Jin and Sun (2004), and Daniel (2003) potential 
energy of a crack body is denoted by 
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cl/r 
cl A 


where 7U is the potential energy of body. The potential energy equals strain energy minus work done. The force induced by 
external agencies is given as 

71 — u — F 


where u = the strain energy stored in the body and 
F = work done by external forces. 


G = J = 


d7T 

~dA 


where G is the rate of change in potential energy with the crack, 
A = crack area, 

u = strain energy stored in the body and 
F = work done by external forces. 


x(l) = JJ v \y /.\ ( / \ ' J* T, tt j d c (1) 

A) c, 

TC = potential energy of cracked solid, 

W — strain energy density 

T j = traction on the boundary segment 
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dA - rdrdO 

\ = A) - ^ 

Substituting equation (6) into equation (5) hence integration along the contour C 0 is reduced to c 


7j = (7 tj n . =0 on the part of the contour C 0 


G = JJ ^-dxdy— J 77 < ^-dc + JJ ^-dxdy— — JJ wdxdy 

. dx dx . dl dl . 

A c A,, A h 

Using the Divergence Theorem and considering the first term of equation (7) 
gives below 

|| —dxdy = | wdy + | wdy 

dx 

c c h 

Substitute equation (8) into equation (7) 

G = | wdy + | wdy — | 77 ^-dc + 1 | dxdy — — 1 | wdxdy 


du 


dw 


d 


c i, 


dx dl 


dl' 


w = B{l)w{X — Z, F) = B{l)w{x, _y) 

h —h 

| wdy = | B(l )vv(— h, y)dy + | b( 1)w(h, y)dy 


JJ ——dxdy — JJ B'(l y )dxdy 


Dl 


Considering the last term in equation (5) 


d_ 

dl 


|| wdxdy 


x = X-l 

w = B(i)w(x - l,Y) = B(l)w{x,y) 


dl 


|| wdxdy = J j m 7- 
A A /->0 A/ 


JJWi + A/ )w(x — A/, v) — B(l )h{x, y))dxdy 


Consider the dummy variable below, 

X* = X- A/ 

h ( h—Al 

|| w(x — Al, y)dxdy — | | w(x*, y)dx* 

A,, —h\ -h—Al 


dy 


(7) 


( 8 ) 


(9) 

( 10 ) 
( 11 ) 

( 12 ) 


(13) 


(14) 
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h f h-M h h-Al \ 



A h -h\-h-Al -h h J 


h ( h-M h \ 



A h -h\-h-M h-M J 

A/ —» 0 


h 

§w(x-Al,y)lxdy=§w(x, y)dxdy+ K w(— h, >’)A/ + w(h, y)Al)dy 

4 4 ~h 

As A/ —> 0 then, 

B(l + Al)= B(l)+ 

Putting equation (13) and equation (18) into equation (13) gives 

7 h 

— jj Wdxdy = B'(l )JJ w(x, y)dxdy+ /?(/)J*(w'(— h, y)— w'(h, y))dy 

A/, A h —h 

Substituting equation (11), equation (12) and equation (19) into equation (9) 

G = | wdy -1 T j ds = J 

dx 


(15) 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


This is the energetic path independent contour integral near crack./ = 0 for any close contour given that Ci and C 2 are two 
arbitrary contours around the crack tip connected by segment along crack face C 3 and C/ then a closed contour is formed.The 
total J along the contour is equal to the sum of contribution from each segment: 

y = j + ! + j + ! =0 

Cj c 2 c 3 c 4 


C 3 = C 4 = 0 and C, = C 2 

Therefore, any anticlockwise path that surrounds a crack will have same value of integral which is path independent. The integral 
is zero over a closed path, Kanninen and Popelar (1985) 


III. GLOBAL ENERGETIC PATH CONTOUR INTEGRAL 

Anderson (2014), Moran and Shih (1987) investigated on the general J - integral where the simplified energy integral is 
modified. However, global representation of J- integral account for dynamic effect, time dependent material behavour including 
the inertial effect. 



V 

F = energy flux 
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J = lim \(w + T) dy - a ii n —— ds 

c =0 J fa 

C 

W — stress work 


w = 


IJ 

\ a u d£ ij 


T = kinetic energy density 


t = L du > du ‘ 

2 ^ dt dt 
d \ d \ 

d {°ij u i)= " 


=p ld +CT '- 


du , 

dX; 


= T + w 


J cT'jU tj in jds = — w + r)v jinjds 

a.. dt .. 


V — volume 

m j = outward normal to surface dv and 


( 21 ) 


( 22 ) 


V j — instantaneous velocity of dv 


In a case of two dimensional flawed solid materials where the partial separation is along x- axis and the origin at crack tip. Define 
a contour to fixed in space containing advancing crack and bounded by the area A while the crack tip is enclosed by a small 
contour c fixed in size and move with the crack. Therefore, the energy balance law is given as: 


J (T ij u [ m j dc = — J J (w + T)dA — J J [(w + T)vd lJ + cj-u-] m .dc 

a c 

J" <7-UfMjdc = the rate at which energy is input into the body 


d_ 

dt 


J (w + T)dA = the rate of increase in energy of the body 


The last term is the rate at which energy is lost from body due to flux through c 


n. = -nij on C 


F ( c ) = J [o + T ) v $u + &,jUi ] n j dc 


(23) 


The flux does not depend on shape of c: therefore, flux to crack tip is given as F = lim I [( W + T)V §■■ + CT ij U i \n :dc 

c —>0 J J J J 

c 

(24) 

As time increases by dt, fissure extends by dl = Vdt and energy expended 
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is fdt therefore, energy release rate is given as 


J = 


F 

v 


(25) 


substituting equation (24) into equation (25) give rise to a global expression for J-integral 


dll- dll- 
u, = -v— L + —- 
dx dt 


(26) 


Hence, in steady state situation, the second term in equation (26) vanishes and the displacement at fixed distance from 
propagating crack tip remains constant, but close to crack tip displacement changes rapidly with positions (at a fixed time). The 
first term is dominative in all cases in equation (26) 


J - lim 

c —>0 J 




rijdc 


Lim 

c-» 0 J 


(w + T ) dy - a ;j n 


ij j 


du i 

dx 


rijdc 


(27) 


The above equation applies to all categories of material response: elastic, plastic, visco-plastic viscoelastic behavior. It was 
derived from globalize energy balance. 


Equivalent Domain Global Energy Integral 

Reza (2014) and Nguyen (2014) also supported the generalized integral 

(w + r)£..-cr.. dU: 


J = lim f 
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dx 


n j dc 
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2 dt dt 
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J =\ 


dU i X 
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2 j 


du ( 

dx x 


qdc 


(28) 


q = 


i a - 

A V 

du t 
dx j 

1 dx, 


A/ dx. 


jd Xi v 




-+- 


ij dx , ox 


50 _ du, 

+ a<j u - F. —- 

dx j dx x 


r du, 

qdA - cr, —-qdc (29) 

+ J - dx. 

c +c I 


0 = thermal effect, 
F j = bodyforce 
p = plasticityeffect 
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The above general representation of J- integralapplies to all 
categories of material response: elastic, plastic, visco-plastic 
viscoelastic behavior.The global integral permits accurate 
analysis possible for all deformation fields including very 
close proximity of crack tip. 

IY. CONCLUSION 

Elastic plastic fracture dynamics that can be used in material 
design and evaluating the integrity of structurenotably among 
which is J- integral, a parameter that deals with energy 
released rate while stress intensity factor is a local parameter 
that deals with displacement and stress field in the vicinity of 
a fissure. The integral quantifies the magnitude of plastic 
zone effect that affect actual fatigue crack progression. 
Therefore, elastic plastic damage dynamics model for crack 
progression becomes necessary. This gives rise to energy 
flux integral that denotes path independent energy line 
integral and the strength of the singular stress and strain field 
near the crack tip that is a measure of energy available for an 
increment of crack extension or the energy release rate. It is 
also referred as crack driving and extension force. This is 
computed from path independent integral in the 
neighbourhood of flaw tip which has the potential of the 
release of energy from the system per unit area extension of 
crack growth. J integral is useful in numerical calculation of 
stress intensity factors. The global energy line integral 
account for dynamic effect, time dependent material 
behavior including the inertial effect. It is applied to all 
categories of material response: elastic plastic, visco-plastic 
material for a better evaluation. The global integral permits 
accurate analysis possible for all deformation fields including 
very close proximity of crack tip. Fatigue crack growth 
analysis are developed to support economic fail safe and 
damage tolerance philosophies and the execution, 
effectuation of damage tolerance control procedures for the 
workings of structural systems. 
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